Correlated flux noise and decoherence in two inductively coupled flux qubits 
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We have studied decoherence in a system where two Josephson-j unction flux qubits share a part 
of their superconducting loops and are inductively coupled. By tuning the flux bias condition, 
we control the sensitivities of the energy levels to flux noises in each qubit. The dephasing rate 
of the first excited state is enhanced or suppressed depending on the amplitudes and the signs of 
the sensitivities. We have quantified the 1// flux noises and their correlations and found that the 
dominant contribution is by local fluctuations. 
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The presence of low-frequency flux noise in supercon- 
ducting devices has been known for decades. It was dis- 
covered as an excess noise in SQUIDs 1,2 and has recently 
been intensively studied as a source of dephasing in var- 
ious types of superconducting qubits, such as charge,- 
fluxj 4 ^— and phase qubits<2r— The noise spectrum typi- 
cally follows 1/f frequency dependence with a spectral 
density of about 1-40 n^o/Rz 1 ^ 2 at 1 Hz, where $ is 
the superconducting flux quantum. The noise spectrum 
depends only weakly on geometry — samples with a loop 
area of 1 /1m 2 up to a few mm 2 are within the above 
range — and does not show a clear dependence on the 
material used. However, the microscopic origin of the 
noise has been elusive so far. It is crucial to identify and 
eliminate the source of such noise in order to improve 
the performance of these devices; i.e., the sensitivity of 
SQUIDs and coherence of qubits. 

A few recent experimental observations have suggested 
as the source of the flux noise a high density of electron 
spins existing on the surface of superconducting elec- 
trodes. Rogachev et al. attributed the magnetic-field- 
induced enhancement of the critical current in supercon- 
ducting nanowires to suppression of spin- flip scattering at 
the surface,— while Sendelbach et al. directly measured 
the magnetization of SQUIDs at low temperatures as 
well as correlated inductance fluctuations! 13 ' 14 Electron 
spin resonance studies on Si/SiC>2 interfaces^ as well as 
scanning-SQUID measurements on Au surfaces^ have 
also indicated the presence of electron spins. As the mi- 
croscopic origins of these surface spins, theoretical studies 
have proposed localized electrons at disordered interfaces 
between surface oxides and metals/semiconductor*^— 
or at defects in the surface oxides^ 

In the present study, we approached this issue by 
means of dephasing measurements in coupled flux qubits. 
Decoherence in coupled qubits depends on the correla- 
tions between noises in each qubit £lr— It can therefore 
be used to characterize the noise correlations without the 
need for measuring the correlations directly. In the echo 
decay signal of the qubits, we observed the contributions 
of pure dephasing due to 1/f flux noises and evaluated 
the correlations. The results indicated local rather than 



global flux fluctuations in accordance with the surface 
spin model. 

The experiments were carried out using a sample fab- 
ricated by electron-beam lithography and shadow evap- 
oration of Al films, with a 20 nm-thick first layer and 30 
nm-thick second layer, on a non-doped Si wafer with a 
300 nm-thick SiC>2 layer (Fig. [1]). The qubits formed a 
small superconducting loop intersected by four Josephson 
junctions, among which one was smaller than the others 
by a nominal factor of 0.55. The sample consisted of two 
flux qubits, ql and q2, coupled with each other via kinetic 
inductance of the shared part. The area ratio of ql and 
q2 was designed to be 1:3 so that the two qubits could be 
closely biased to their half-integer flux-bias points simul- 
taneously by a global magnetic field. The conditions were 
$cxi/$o = 0.5 and $ e x2/$o = 1-5, where $ eXj is the ex- 
ternally applied flux through qubit j. An additional on- 
chip local flux bias line made of Al allowed independent 
control of the flux biases in each qubit. 

The effective Hamiltonian of a system with two induc- 
tively coupled flux qubits can be expressed in terms of 
persistent current basis as 
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where a X j and a Z j are Pauli matrices, Aj is the tun- 
nel splitting between the two states with opposite direc- 
tions of persistent current along the loop, £j = 2I p j<&onj,j 
is the energy bias between the two states, and I p j is 
the persistent current along the qubit loop. The qubits 
are coupled with a coupling energy J12 = M\%Ip\Ipi, 
where M12 is the mutual inductance. We define the 
normalized magnetic flux in each qubit loop n^j as 
{(^gxj — 0.5<f>o) m °d ^ol/^o- I n the case of an isolated 
qubit, Eqi = Aj and dEoi/dn^j = at n$j = 0, where 
Eqi is the eigenenergy of the first excited state. This is 
the optimal flux bias condition where dephasing due to 
fluctuations of n^j is minimal. 

In our experimental setup, the energies of the first and 
second excited states, Eqi and -E02, in the two-qubit sys- 
tem are measured by spectroscopy. A 5 /is microwave 
pulse is applied to the system, followed by a bias current 
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FIG. 1: (color online) Scanning electron micrograph of the 
sample and a sketch of the measurement setup. Two coupled 
qubits (ql and q2) share a part of their loops with a read- 
out SQUID shunted by an on-chip capacitor. A current bias 
line (lb) and a voltage measurement line (V m ) consisting of 
lossy shielded cables are connected to the SQUID via on-chip 
resistors. Microwave current pulses are fed through an on- 
chip control line (Imw) inductively coupled to the qubit. A 
local flux control line (Idc) is also inductively coupled to the 
qubit and is connected to a battery-powered current source. 
The sample chip is enclosed in a copper-shielded box, and 
all the wires are electrically shielded as well. CPF and RCF 
denote a copper-powder filter and RC low-pass filter, respec- 
tively. Global magnetic flux bias is applied with an external 
superconducting coil connected to a battery-powered current 
source. The sample is cooled to 20 mK in a dilution refriger- 
ator magnetically shielded with three /U-metal layers at room 
temperature. Inset: Schematic of two coupled qubits. This 
system consists of three segments: a part of ql (LJ1), a part 
of q2 (LJ2) and a common part (L12). 



pulse of the readout SQUID. When the microwave fre- 
quency hits a transition of the system, the excitation is 
detected as a change in the SQUID switching probabil- 
ity P sw i2£ Figure [2] shows the results as a function of the 
global flux bias. By sweeping the global flux bias, n§\ and 
n^2 vary according to the equation ri0i = 0. 3371^2 + n^s, 
where n^g is an offset between the optimal flux bias 
points of the two qubits. For n^s — 0, the optimal bias 
points of ql and q2 are aligned, while the optimal point 
of q2 shifts toward the left with decreasing n^s ■ The en- 
ergy levels show anticrossings due to the strong inductive 
coupling between the qubits. 

The data fit well with the calculated eigenenergies, in- 
dicating that the system can be simply described with 
two qubits and a fixed coupling between them. The fit- 
ting parameters obtained from the least-squares method 
are as follows: A x /h = 6.56 ± 0.02 GHz, A 2 /h = 
5.29 ± 0.03 GHz, I pl = 125 ± 1 nA, I p2 = 136 ± 1 nA, and 
Jx-z/h = -1.20 ±0.02 GHz. By partially differentiating 
-Soi (7i0i j 7102) with respect to n^i and ri02, we also obtain 
energy sensitivities of the first excited state in the two- 



FIG. 2: (color online) Spectrum of the two-qubit system as a 
function of the global flux bias. Panels (a)-(c) are for differ- 
ent flux bias offsets n^s- The open circles are measured data 
corresponding to £01 and E02, the energies of the first and 
second excited states. The red solid curves represent calcu- 
lated energy levels using the two-qubit Hamiltonian. 
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FIG. 3: (color online) Sensitivities of eigenenergies E01 to flux 
noise in each qubit loop as a function of the global flux bias. 
Panels (a)-(c) are for different flux bias offsets n$s- The black 
solid curves are for dEoi/dn^i, and the red dotted curves for 
dE i/dn$2- 



qubit system to flux noises in ql and q2; i.e., dEoi/dn^i 
and dEoi/dn^ (Fig-El), respectively. 

Dephasing of the first excited state is character- 
ized by Hahn echo measurement with a sequence of 
(7r/2)-(7r)-(7r/2) pulses.— The echo decay curve P sw (t) 
contains both relaxation component exp(— Fit/2) and 
pure dephasing component {e lv ^)E such that P sw {t) oc 
cxp(— T it / 2) {e 1 ^}) e ■ The first term is determined inde- 
pendently by energy relaxation measurement, while the 
latter is fitted well by a Gaussian decay, exp{ — (T^ E t) 2 }<^ 
Figure [4] shows T^ E obtained for three different n$s- 
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FIG. 4: (color online) Pure dephasing rates (circles) deter- 
mined by spin-echo measurements as a function of the global 
flux bias. Panels (a)-(c) are for different flux bias offsets n^s- 
The red solid curves represent a fitting considering the ampli- 
tudes of the flux noises in each qubit and their correlations. 
The blue dotted curves were calculated by subtracting the 
contribution of the correlations from the red solid curves. 



The sample has relatively good coherence: the smallest 
dephasing rate T 9 E = 0.153 x 10 6 s" 1 is observed for 
n<j,s = —8.75 x 10~ 3 and at n^\ — —2.85 x 10 -3 , where 
T 1 = 0.148 x 10 6 s -1 . The decay of the echo signal by a 
factor of 1/e takes place in T 2oc ho = 4.4 /xs. 

We have also examined the SQUID-bias-current de- 
pendence of (e^V^E and confirmed that the noise from 
the bias current line does not significantly contribute to 
the pure dephasing.- Because of the strong n^j depen- 
dence of r^ E , we can also rule out possible contributions 
of charge and critical-current noise to the dephasing. 

The observed Gaussian decay implies the presence of 
low-frequency energy fluctuations with a 1/ / spectrum, 
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It also requires that the high-frequency tail extends up to 
the angular frequency range comparable to T 9 E , which 
is ~ 2ir x 10 MHz in the present sample. The echo de- 
phasing rate is then expressed as T 9 E = y/ Ae 01 ln2/?i. 28 
As the origin of the energy fluctuations, it is most nat- 
ural to assume the presence of 1// flux noises in each 
qubit: S n ^(u) = A„ w /M (j=l,2)A Here we also con- 
sider cross correlations between the two flux fluctuations, 
Sn^,i(t) and 8n^%{€). Following the linear cross approxi- 
mation discussed in Ref. 24, 
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This form of cross-correlation spectrum is justified when 
the 1// noise is produced by many two-level systems and 
their coupling strengths to the qubits and switch rates are 
uncorrelated. For 1/f charge noise, such cross correla- 
tions in charge noise have been observed using two single- 
electron transistors.— Now we fit the observed Gaussian 
decay rate T\ g E using 
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(red solid curves in Fig. U]) • All the data points of T 9 E 
for different flux biases arc fitted with a single set of pa- 
rameters, A n<pi = [(2.32 ±0.01) x 10" 6 ] 2 , A n4>2 = [(2.76 ± 
0.01) x 10- 6 ] 2 , and A n ^ 2 = [(1.49 ± 0.01) x HT 6 ] 2 , 
which confirms the validity of our assumption of dephas- 
ing induced by 1// flux noises. Moreover, the noise am- 
plitudes, similar to the previously reported values, are 
determined with high accuracy. We also note that the 
T 9 E data cannot be fitted well if A n4>in4i2 is set to zero. 

The small ratios A n ^jA n ^ = 1.42 and 
An^n^ / A nipl = 0.42 are clear proofs that global 
flux noise is not a dominant source for the dephasing 
in the first excited state: global flux noise would give 
rise to A n4>2 /A n<t>1 = 9 and A n4>in<f>2 / A n<t>1 = 3 due to the 
difference between the areas of ql and q2. 

On the other hand, we obtain a non-negligible amount 
of An, in , 2 , which indicates noticeable correlations be- 
tween Sricjixif) and Sn^if). In Fig. 01 dephasing rates 
calculated without taking into account the correlation 
term arc plotted as blue dotted curves. The observed 
dephasing rate is enhanced (suppressed) when sensitiv- 
ities BEgi/dn^i and dEoi/dn^ have the same (oppo- 
site) signs. For example, at around n^i = — 6 x 10~ 3 for 
n<ps — —8.75 x 10~ 3 , reduction of the dephasing rate due 
to the correlation is observed. 

To qualitatively account for the flux noise amplitudes 
A n , A n , and especially the correlation term A n n , 
we introduce a simple model where flux noises are gen- 
erated by a number of microscopic sources, described as 
fluctuating magnetic dipoles, scattered over the sample. 
The two-qubit system can be divided into three segments 
(see the inset of Fig. [IJ: a part of the ql loop and its 
junctions (LJ1), a part of the q2 loop and its junctions 
(LJ2), and the common part of the two loops (L12). Cou- 
pling between each dipole and a qubit reaches a maxi- 
mum when the dipole sits on the surface of the super- 
conducting loops, while dipoles away from the surface do 
not couple effectively^ For such sources on the surface, 
we further apply the following approximation: Sources 
along LJ1 (LJ2) couple exclusively to ql (q2) and those 
along L12 couple equally to both ql and q2. For the 
latter, Sn^i and bn^i have full correlations. 

Summing up all the contributions of dephasing from 



LJ1, LJ2, and L12, Ae 01 can be rewritten as 
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The noise amplitudes in each part of the system are cal- 
culated as A™ 1 = (1.77 x 1CT 6 ) 2 , A^f = (2.32 x 1(T 6 ) 2 , 
and A^ 12 = (1.49 x 10~ 6 ) 2 , which represent the distribu- 
tion of flux noise sources on the loops. Two points are 
worth mentioning here: (i) There are two parts of the 
system, (L12 + LJ1) and LJ2, each consisting of four 
Josephson junctions and a similar length of supercon- 
ducting loops. The noise amplitudes in the (L12 + LJ1) 
part are calculated as A^ 1 + A^ 12 = (2.32 x 1CT 6 ) 2 , 

which coincides with A^f 2 . This result is consistent with 
the model and the assumption of the local flux noises, 
(ii) The third term in Eq.([S]) originating from L12 can 



be cancelled when the sensitivities satisfy the equation, 
dEoi/dn^i = -dE Q1 / Orig- 
in conclusion, we have studied dephasing in two induc- 
tively coupled flux qubits. The dominant source of the 
dephasing is found to be low-frequency 1// flux noises. 
At the same time, the local, rather than global, nature 
of the flux noises is revealed in the correlations between 
the noises in each qubit. The results agree with a model 
in which flux noise sources are distributed on the surface 
of superconducting loops. 
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